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Abstract 

A new expansion is established for the Green's function of the elec- 
tromagnetic field in a medium with arbitrary e and /i. The obtained 
Born series are shown to consist of two types of interactions - the 
usual terms (denoted V) that appear in the Lifshitz theory combined 
with a new kind of terms (which we denote by Q) associated with 
the changes in the permeability of the medium. Within this frame- 
work the case of uniform velocity of light (e/i = const) is studied. 
We obtain expressions for the Casimir energy density and the first 
non-vanishing contribution is manipulated to a simplified form. For 
(arbitrary) spherically symmetric we obtain a simple expression for 
the electromagnetic energy density, and as an example we obtain from 
it the Casimir energy of a dielectric-diamagnetic ball. It seems that 
the technique presented can be applied to a variety of problems di- 
rectly, without expanding the eigenmodes of the problem and using 
boundary condition considerations. 



1 Introduction 

A theory describing the fluctuations of electromagnetic fields in dielectric 
medium was developed long ago by Lifshitz, Pitaevskii and others, and is 
now classical textbook material ^ [3[]. Among other things the theory 
permits efficient and general treatment of the Van der Waals interaction and 
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calculation of Casimir energies. This theory was developed mainly for cases 
where the medium has arbitrary unhomogenous permittivity e(r, uj) while 
the material considered is not magnetoactive (with uniform jj). In another 
classical work, Balian and Duplantier |4[] derived an expansion for the Green's 
function of the electromagnetic field in vacuum, but in the presence of a 
conducting boundary surface. The terms in this expansion are related to 
the geometry of the surface, and are given as explicit integrals, which differ 
from terms encountered in |]J . This expansion was then applied to the study 
of the statistical properties of the electromagnetic eigenmodes and led to 
an expansion for the Casimir energy under arbitrary conducting boundary 
conditions at zero and finite temperatures ||. 

Although having no direct relation to the dielectric material problem, for 
the special case where e/x = const (this condition will be referred to as the 
uniform velocity of light - UVL in short) there is such a connection. This 
was first apparent in the spherical problem where the UVL scenario was 
first introduced in the present context by Brevik and Kolbendsvent |§. The 
problem of a dielectric-diamagnetic ball in a medium under the UVL condi- 
tion is evaluated by expanding either the Green's function or the appropriate 
characteristic functions for the modes in terms of the diluteness parameter 
£ = (n — \i ') I \\i + //) (where /i and y! are the permeabilities inside and 
outside the ball respectively). It was shown that one can calculate the term 
proportional to £ 2 || exactly by means of an integral over the free Green's 
function of the Helmholtz equation. While introduced as a formal trick, it 
coincides exactly with the two scattering result in || up to multiplication by 
£ 2 . This leads us to develop a general expansion for the Green's function in 
dielectric-diamagnetic media. 

The expansion that will be presented coincides with the approach of Lif- 
shitz when /i is constant, and bears strong formal similarity to the multi- 
ple scattering expansion when considering dielectric bodies immersed in a 
medium in which UVL holds. It is interesting to point out that while the 
parameter encountered in many calculations in the literature of the Green's 
function in the latter case is £, the natural parameter appearing in the present 
expansion is log[(l — £)/(l + £)]. 

We start by deriving a Born type series for the Green's function in a 
medium with arbitrary e and \x. The series are shown to consist of two types 
of interactions - the usual terms (which will be denoted V) that appear in 
the Lifshitz theory combined with a new kind of terms (which we denote 
by Q) associated with the changes in the permeability of the medium. The 
V terms are proportional to (1 — e/i) and thus, when UVL holds only Q 
terms contribute to the Green's function. Next we investigate the properties 
of the Q terms are investigated and it is shown that in the UVL case the 
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magnetic and electric Green's functions are closely related - a feature which 
is not present in general settings of dielectric medium. We proceed to use the 
Green's function for the calculation of the Casimir energy of the medium. We 
study in detail the contribution of the second term in the expansion to the 
energy density of the electromagnetic field, and obtain an integral formula. 
For the case of an arbitrary radially dependent /i this contribution can be 
considerably simplified. As an example, we use the radial formula to calculate 
the Casimir energy of a dielectric-diamagnetic ball and show that the exact 
energy to order £ 2 is easily recovered. 



2 Born type expansion 

The statistical properties of the electromagnetic field in a medium are de- 
scribed by the appropriate photonic Green's function. The electromagnetic 
fields are derived from the electromagnetic potentials A a , a = 0, ..,3. (It 
is convenient to work in the gauge A = 0.) The Green's function is now 
defined by: 

V ik (X 1 ,X 2 ) = % < TA i (X 1 )A k (X 2 ) >, (1) 

where Xi,X 2 are 4-vectors (X®, --Xf) and i, k — 1, ..3. The angular brackets 
denote averaging with respect to the Gibbs distribution and T denotes the 
time ordering operator. In general we will be interested in the retarded 
Green's function, D ik , which is defined by: 

D (X X ) = i 1 < MXi)A k (X 2 ) - A k (X 2 )A i (X 1 ) > t x <t 2 
ik{ - u 2 ' l0 otherwise ^ ' 

This function, when Fourier transformed with respect to the time difference 
is actually the generalized susceptibility of the system [jjj |j| . It is known that 
in a medium with a given permittivity tensor and permeability tensor fj,^, 
D satisfies the equation: 

U) 2 

V x lTV x D + — eD = -4nhl5(r - r') (3) 

c 2 

where I is the 3-dimensional unit matrix. In what follows we shall work in 
units where c = h = 1. 

It will be convenient to introduce the following notation, which is suitable 
for some of the calculations ahead: Let a be a 3-vector; we define the cross 
product operator with a as [a]i k = €ij k aj. This operator acts on vector valued 
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functions byQ [a] f = ~~ct x / . Thus, for example, the curl operator is denoted 
[V]. Furthermore we adopt the convention that the same notation is used 
for the operators and their kernels. With this notation the equation for 
the Green's function (multiplied by 4tt) of dielectric media has the following 
form: 

([V]^ 1 [V]-uj 2 e)D = -I (4) 

In the following we assume that /x(r) is a scalar function (of course, in the 
general case both /i and e are tensors). 

To write this equation in a manner convenient for the following we multi- 
ply by /i and use the fact that /i[V]/i _1 = — [Vlog/i] + [V]. Thus obtaining: 

([V] 2 - [Vlog/x][V] - u?efi)D = -//I (5) 

The inverse of the operator ([V] 2 — uo 2 ) is well known and will be denoted 
D . The kernel of D , which is the Green's function in vacuum, is given by 
the formula: 

D (oo; r, r') = - (l+ ® VJ g (r - r'), (6) 

where go is the Green's function for the Helmholtz equation: 

Ag Q + u 2 g Q = -5(r-r'), (7) 

given by: 

g (r-r')= 1 exp^l^l). (8) 
47r|r — r'| 

Now we wish to use the known D$ to express D via a Born type series. 
Denote P = V + Q where: 

V = uj 2 1{\ - lit) (9) 

and 



Q = -[Vlog(/x)][V] (10) 

2 The convenience of this notation will arise from the necessity to evaluate expressions 
such as Tr(a x b x cx). Written in this form one has to keep track of the order of the 
multiplication because of the nonassociative nature of vector products. Working simply 
with matrices such as [a] will prove to be very transparent. 
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Then: 



D = ((/ - PD )D - 1 )- l fi = D (I - PD )- l fi (11) 

Thus we get the following formal series for D: 

D = D of i + D (V + Q)D ot i + D (V + Q)D (V + Q)D of j, + .... (12) 

Thus the series for D consist of two types of interaction. The V terms 
are well known and can be found in M. However, to our knowledge the 
magnetoactive terms Q are new in this context and are related to the changes 
in the permeability of the medium. 



3 The case of uniform velocity of light 

Uniform velocity of light in the medium implies efi = I. This considerably 
simplifies the expansion fll2"|), since it just leaves the Q scattering terms. We 
are then left with the expansion: 

D — DofM + DqQDqh + DoQDoQDofM + .... (13) 

In the following we are going to apply Q to the free Green's function D Q . For 
this purpose it is convenient to introduce?] 

G = [V]A, = [V<7o] (14) 

Thus the action of Q on Dq is given by 

QD = -[V log n)G (15) 



So that ([13]) can be written in the form: 

D = D of i- DofVlog^Go/i + ...(-l) n - 1 J D ([Vlog yU ]G )> - +.... (16) 

The terms in (O) can be written as explicit integrals using relation (|A.5|) : 

D(r,r') =/i(r') J D -Mr , )/^o(r,r 1 )( 7o (r 1 -r')[Vlog/i(r 1 )][u ri ,,]+ (17) 
t*(?) J D o(?,*i)9o(ri - r 2)^ ( r 2 - r')[Vlog/i(r 1 )][u rir2 ][Vlog/i(r 2 )][u r2r /] - .... 

Convergence of the series relates of course on the nature of /i(r). For specific 
cases, such as of a body immersed in dielectric-diamagnetic medium the 
expansion will have a natural parameter, as will be explained in section 4. 

taking the rotor of Dq simplifies this function considerably since we are only left with 
first order derivatives of the scalar Helmholtz propagator: indeed Cijidjdidk = 
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3.1 The magnetic Green's function and its relation to 
the electric Green's function 

Another quantity of interest is the magnetic Green's function, given by: 
x D x V . In this section we show that the expansion for the mag- 
netic Green's function is similar to that of the electric Green's function, the 
main difference being alternating signs of the terms. This property is impor- 
tant for calculation of the Casimir energy, since it cancels odd terms in the 
expansion. 

Indeed let us write 

D = (/ + D [Vlog/i][V])- 1 IV = A)(/+ [Vlog/zJCorV (18) 

Thus, 

V x D x V = (19) 
[V]([V]Z^ = [V]([V]^ (/ + (D [V log f,][V]Y)- l Y = 
[V](/ + J Do[Vlog/i][V])- 1 (- J D [Vlog /U ]/i+ [V]£M = 
(I + G [ V log n] ) - 1 ( -Go [ V log n] n + u 2 D 0fi - fj) = 
^(J + GofVlog/i])- 1 ^-^ 

(Here we used the elementary property (1 + AB)~ 1 A = A(l + BA) -1 .) We 
can now write the expansion of the magnetic Green's function in the form: 

X? x D x V = (20) 
-y. + u 2 D y + u 2 {D [V log fjtjGoYfi + ...u 2 {D {[V log^]G ) n )V + 

Comparing this expansion with the one for D (Eq fllBD ), we see that the 
magnetic Green's function has the same type of expansion as D f , but with 
constant signs, owing to [V log fi] being the only antisymmetric operator in 
the expansion. A similar symmetry was obtained in a different way for the 
Green's function in vacuum with conducting boundary conditions in [[| and 
its implication on the distribution of eigenmodes is discussed. The impli- 
cation of this symmetry to the calculation of the Casimir energy, is that in 
the series for the energy the main term one has to consider is the second 
order term in the expansion. This follows since the zeroth term is usually 
subtracted as the free problem solution while the first term cancels between 
the electric contribution and the magnetic contribution. 

3.2 The symmetry when relaxing the UVL condition 

What can be said when keeping, in addition, the V terms? Can we still 
have cancellation or similarity between the electric and magnetic Green's 
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function? Cancellation in the spirit of the preceding section can be achieved 
by expanding as follows: Let D e be the solution of 

{[V][V]-oo 2 ey)D e = -l (21) 

And accordingly G e = [V]-D e . Then it is possible to expand D in terms of 
D e . Using the fact that D e is a solution for a non magnetoactive problem and 
therefore has the symmetry property D^r, r') = Djy(r', r). One may show 
that a there is a again a symmetry between the expansions of the magnetic 
and electric Green's functions, similarly to the one discussed in the previous 
section. 



4 Dielectric-diamagnetic body immersed in a 
medium 

In this section we consider a body B with constant dielectric properties e, fi 
immersed in a medium with different dielectric properties e',//. Let <9B 
denote the boundary of the body, and let s be coordinates on this boundary. 
In this case the Q tensor, as a function of its first variable, will be supported 
on the boundary of the body (more generally, on <9B x R 3 ) and has the 
following structure: 

Q(s,r) = log(4)[n(s)]G (s,r) (22) 

where n(s) is the unit normal to <9B at s. 
When UVL holds, the expansion for D is: 

D(r, r') = Mr') A) - log(£) J m d Sl D (r, s 1 )[n(s 1 )]G (s 1 , r') + (23) 
log(^) 2 / 9B c/sirfs 2 D (r,s 1 )[n(si)]G'o(s 1 ,S2)[n(s 2 )]G'o(s2,r / ) - .... 

This can be further simplified using (|T7|) : 

D(r, r') = Mr') A) - log(f )/i(r') f m d Sl D (r, Sl )g' (s - r') [n s ] [u sr ,] + (24) 
M(r') log(^) 2 / aB dsids 2 D (r,Si)^(si - s 2 )^(s 2 - r')[n Sl ][u Slfl2 ][n S2 ][u S2r /] - .... 

Thus the Green's function is expanded in terms of the parameter: log(^r). 
It is interesting to relate this parameter to the parameter £, which appears 
in the literature when considering the case of dielectric-diamagnetic ball 
or cylinder. To see the relation we expand: 

log ^ = log = -2£ - 2^ + ... (25) 

u 1 + 1 3 
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Hence the two expansions can now be related by assuming £ small. In par- 
ticular the second term, proportional to (log^-) 2 is exactly one forth of the 
coefficient of £ 2 in the common expansions of the Casimir energy 

5 Energy density from Green's function 

Fluctuation dissipation theory says that (with our definition of D), at tem- 
perature T = 4 the correlation functions of the fields Ai are related to the 
generalized susceptibility D by 0: 

(A(r)^(r')). = -z47rcoth(^)(A,-(r,r') -%(r',r)) (26) 
Now one can check that 

(EiE k ) u = u^AiA^ (27) 

and that 

(BiB^ = V x (AA k ) w x V. (28) 

Note that there is no minus sign before the correlator of the electric fields, 
this can be verified by returning to the time parameter and using E = —A. 
Thus the energy density per du of the electromagnetic field at temperature 
T is given by 

p T (r,u,) = -coth(^)p(r,u;) (29) 

where 

p(r, uj) = -^Limit r /^ r ImTr[cj 2 e(r)D(r, r') + xDxV] (30) 

4:71 M 1 ") 

This expression for p is usually divergent and should be regularized. The 
natural regularization scheme in this problem would be to subtract the energy 
density of the spatially homogeneous problem for this result. 



6 Expansion for energy density and density 
of modes in a medium with UVL 

As was previously shown, in the UVL case there is a symmetry between 
the electric and magnetic Green's functions. As a result the calculations are 
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considerably simplified. In the expression for the energy density the zeroth 
order terms are cancelled by the vacuum D function, and moreover we are 
left solely with the even terms of the expansion. All the higher order terms 
are given by explicit integrals in which the limit r' — ► r can be taken at 
the outset. Thus, for example, the first term to be considered in the energy 
density is: 

p T ( 2 )(r,^) = -coth(^)p( 2 )(r,c) (31) 

Where 

p( 2 )(r,^) = ^uHmjg^ - r 2 )g' (r 2 - r) (32) 
xTr([Vlogp(r 1 )][u rir2 ][Vlogp(r 2 )]G (r 2 ,r) J D (r,r 1 )) 

To find the density of the Casimir energy with respect to uj this expression 
is to be integrated over the volume considered. It is convenient to integrate 
first over r (see 0]). In our notation: 

J drG (r 2 ,r) J D (r,r 1 ) = [V r2 ](/+ ^V ri ® V ri ) J dr g (r 2 - r)g (r - n) 

(33) 

Using the identity: 

J drs-o(r2 - r)g {r - n) = ^d w g G {T 2 - ri) (34) 

We obtain the result: 

J dr G (r 2 , r) J D (r, n) = ^d u G (r 2 - n) (35) 

Namely we have the operator identity, Go^o — Y~d u Go. This can be done at 
any order in the perturbation series so that we obtain the following general 
expression: 

p( 2 ™V) = ^ 2 Im/Tr([Vlogp( ri )][u rir2 ] • ■ • [V logp(r 2n )][u r2nri ]) (36) 
x#o( r i - r 2) ■ • ■ £o( r 2n-i - r 2n )^d^o(r 2n - n) 

In principle the Casimir energy can be obtained from these expression to any 
order and for any p. However, in most cases this is impossible to carry out, 
and one has to be content with a first few orders. 
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7 Casimir energy: second order results 

In this section we focus on the second term in the expansion for the density 
(|36f) . We will obtain a simplified expression for this term, however, one must 
bear in mind that for any specific problem it is necessary to check whether 
this term is indeed larger then the other terms in the expansion. To cast the 
second term in a more convenient form, it is useful to introduce the function: 

I g2io;|r2-ri| ^ 

92 ( ri - r 2 ) = g' Q {v x - r 2 )-cU(r 2 - r x ) = ^^j (^r^| " ™)> 

(37) 

so that the second order contribution to the energy is given by fl36|): 

P (2 V) = ^ 2 Im / TlCtVlog^rOl^ltVlog^ra)]^]) (/ 2 (n - r 2 )(38) 

Using the relation ( A . 1 (J| ) 

7Y([a][b][c][b]) = 2(a-b)(b-c) (39) 

We obtain 

p {2) (uj) = ^uj 2 lm y'(Vlog / u(r 1 ) • u rir2 )(Vlog/i(r 2 ) • u r2ri )5f 2 (ri - r 2 ) (40) 

This is our final result for the first contribution to the Casimir energy for 
arbitrary /i. For simple geometries this can be still simplified. For example, 
if we assume that fi is a function of one coordinate only, namely /i = fi(z), 
Then [Vlog/x] = (log/x(2;))'[z] and the second term is 

PT (2) M = -i^ 2 coth(^)Im J 2(z- u rir2 ) 2 ^ 2 (r 2 -r 1 )(log / i(^ 1 ))'(log/i(z 2 ))'(41) 

The energy density per unit area for this configuration can be now obtained 
by integrating this expression over r 2 and the z coordinate of rj.. 



7.1 Second order term for spherically symmetric /i 

For cases with high symmetry the equation (|40D for the energy density can be 
further simplified. Indeed, let us assume that we are dealing with a medium 
with spherical symmetry. This type of mediums can be encountered in a 
variety of problems from QCD bag models to cosmology. In our formalism 
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this amounts to consider a magnetic permeability \x which is a radial function. 
Thus [Vlog/i] = (log/i)'[u 0r ]. The second term ( fiO|) is then given by 

P (2) M = -£lm/Tr([u 0ri ][u rir2 ][u 0r2 ][u r2ri ]) (42) 
x^ 2 (ri - r 2 )(log/i(|r 1 |)) , (log y u(|r 2 |)) / 
= ~%^ m - / 2(uo n • u rir , 2 )(u 0r2 • u rir2 ) 
x^ 2 (ri - r 2 )(log / u(|r 1 |))'(log / u(|r 2 |)) / 

To carry the integration further, we change to bipolar coordinates u = |r 2 — | 
and v = |r 2 | and r = |ri|. Thus: 

(v 2 — u 2 — r 2 ) (r 2 — v 2 — u 2 ) (u 4 — (v 2 — r 2 ) 2 ) 
2ru 2vu 4u 2 rv 

(43) 

and the integration becomes: 

Im/ °°47rr 2 dr(logMr))7 2 ^dndy ^^p^g 2 (u) (log ^(v))' (44) 



Where the integration domain is such that u, v and r can form a triangle. 
More specifically, the integrations are given by: 

/ rr j-r+u poo t-u+r \ / / 4 _ / 2 _ r 2\2\ \ 

7(r) = U du y dv + jduj dvj H L- ii^( M )(log / i( V ))'J(45) 

Finally, the second order contribution for the Casimir energy in this case is 
given by the formula: 

p T V)(u) = -47rcu 2 coth(^)Im/ oo dr/(r)(log / u(r)) / (46) 

This equation is our final formula for the Casimir energy of a UVL config- 
uration with arbitrary radially symmetric fi. As an example, in the next 
subsection we implement it for a one-line calculation of the Casimir energy 
of a dielectric diamagnetic ball. 

7.2 Example: The dielectric-diamagnetic ball 

It was first shown by Boyer H, that, contrary to intuition, the Casimir force 
on a conducting spherical shell is repulsive. A large number of subsequent 
calculations by various techniques || ^ verified Boyer's result. In the follow- 
ing, we show that the calculation of the Casimir energy of a ball of radius 
a immersed in a medium, under the UVL, condition is immediate using the 
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technique developed above. In this case \i is of the form: /i(r) = \i if r < a 
and /i(r) = jj! if r > a. Let k = log -4- Then from ([!(]) it is immediate that 



the density of states is of the form 



p T ®(u) = -4W 2 coth(^f)Im(/« 1(a)) 



(47) 



Now explicitly calculating 1(a) we have: 

= Jo ^ 

Performing the integral in fl48f) and substituting in (|47"1), we have: 

1 + e 4iaw - 4ie 4iaw acj - Ae iiaw uo 2 a 2 ' 



(48) 



^(w) = _lm 



8tt 



(49) 



(Where we used ( |25|) to write the result in terms of £.) Taking the imaginary 
part of the expression on the right gives one a simple way of eliminating the 
divergence in this problem, which arises upon integrating over frequencies. 
This simply eliminates the term proportional to 1/87T in the energy density. 
Alternatively, this term is independent of radius and thus can be eliminated 
by subtracting the unconstrained system where a — > oo. In any case we 
recover the result obtained in || by mode summation. 

At zero temperature, the integration over frequencies is conveniently carried 
out on the imaginary axis to yield the known result to order £ 2 : 



E 



f 2 

dcu—e- 4alujl (1 + 4aU + Aa\u\ 2 ) 



2al67r 



(50) 



Since the energy is a positive quantity decreasing with radius the Casimir 
pressure on the boundary will be outward. At finite temperatures, one has 
to take into account the coth(^) factors that appear in (|29|). This is conve- 
niently done by changing the integration to summation over the imaginary 
Matsubara frequencies. In |T(| this was done using the expression (|49|), with 
the result : 



E c (a,T) 



Tj 2 , rg 2 

4 ~+~ 2(e fk " rT -l) 



4 a7r T 2 g 2 e 8 ^ T 



(51) 



2, r 2 T 3£2 



This result is exact (to order £ 2 ) for all temperatures. It is easy to check that 
the pressure on the boundary remains outward at all temperatures. 
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8 Concluding remarks 



In this paper we presented a technique appropriate for evaluation of the 
Green's function of dielectric media, especially in cases where uniform veloc- 
ity of light (UVL) holds. The perturbative technique developed allows one to 
deal with the Casimir effect for geometries not studied in other approaches. 
It was shown that for the UVL case the magnetic and electric Green's func- 
tions are closely related - a feature which is not present in general settings of 
dielectric medium. Due to this property the first contribution to the energy 
density of the electromagnetic field comes from the second term in the expan- 
sion. We have shown that equation (|I0D is an explicit integral representation 
of this contribution. For cases of spherical symmetry we derived a simpli- 
fied formula, namely (pEED, for the energy density of a configuration with an 
arbitrary radial /i. From this formula the energy of a dielectric-diamagnetic 
ball follows easily, without recourse to expansion of the modes themselves or 
boundary condition considerations, other problems can now be also accessed 
in the same manner. 
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Appendix: Identities 

For the convenience of the reader we list here a few simple relations: 
1. 

[b][a] = a®b-a-M (A.l) 
This is easy to see by working with indices: 

([b][a]) im = tijkbjtkimai = (5u5 jm - 8 im 8ji)aibj = a,,b m - 5 im aibi (A.2) 

2. 

G (r,r') = ^(r-r')[ U ^] (A.3) 
where u rr / is a unit vector in the direction from r to r'. Indeed: 

G„(r, r') = [V]A) = -M<7o(r - r')I = -[V (r _ r ^ (r - r')] = ^(r - r')[u rr ,] 

(A.4) 
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3. From these two we conclude that: 

QD = [Vlog/x(r)]G (r,r') = ^(r - r')[V log^(r)][u rr ,] (A.5) 

4. Symmetry properties 

(z) D 0ij (r, r') = D 0ji (r', r) (ii) G 0ij (r, r') = G Qji (J, r) (A.6) 
(Hi) G mj (r,r') = -G Qji (r,r') = -G 0ij (r',r) (iv) [a] = -[a]* 

Remark: Property (i) is true in the form Dij(r, r') = Dji(r',r) for any D 

which is a susceptibility function for non-magnetoactive medium |l[ . 

5. 

A x V = ([V]A7 (A.7) 

6. 

V x D x V = u 2 D Q - I (A.8) 

Indeed: 

^xD xV' = MiMWY = [V]G = [V][V]A> = -I + ^Do (A.9) 

7. Using the relations (a ® b) (c (g> d) = (b • c)a Cg> d and Tr(a ® b) = (a • b) 

we get 

Tr([uo ri ][u rir2 ][u 0r2 ][u rir2 ]) = (A.10) 
Tr((u I)(u !)) = 

2(ll0ri ' UriraJvllOra ' U rir2 ) 
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